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Incoherent pion photoproduction on the deuteron is studied in the first resonance region. The 
unpolarized cross section, the beam asymmetry, and the vector and tensor target asymmetries are 
calculated in the framework of a diagrammatic approach. Pole diagrams and one-loop diagrams 
with NN scattering in the final state are taken into account. An elementary operator for pion 
photoproduction on the nucleon is taken in various on-shell forms and calculated using the SAID 
and MAID multipole analyses. Model dependence of the obtained results is discussed in some detail. 
A comparison with predictions of other works is given. Although a reasonable description of many 
available experimental data on the unpolarized total and differential cross sections and photon 
asymmetry has been achieved, in some cases a significant disagreement between the theory and 
experiment has been found. Invoking known information on the reactions — > ir°d and — > np 
' we predict the total phot oabsorpt ion cross section for deuterium. We find that our values strongly 

^vQ , overestimate experimental data in the vicinity of the A peak. 
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I. INTRODUCTION 



Comprehensive measurements of total and differential cross sections of inclusive, coherent, and incoherent ir° 
J£Z ' photoproduction from the deuteron in the energy region from 140 to 792 MeV were carried out at MAMI 

mi- it 

^ was found that the coherent data are in good agreement with theoretical predictions. However, in the case of the 
incoherent cross sections the situation was found to be much less satisfactory. The theoretical predictions from Refs. 
I> HQ in the A region exceeded the experimental data significantly. 

The above mentioned disagreement may be indicative of shortcomings in the approaches of Refs. HQ. The model 
developed in Ref. Q seems to be oversimplified because it takes into account the pole diagrams only. It is known that 
nuclcon-nucleon final state interaction (FSI) is extremely important in incoherent pion photoproduction especially for 
small pion ang les (see, e.g., Refs. HUH). Alt ho ugh FSI was incorporated in the model of Ref. H > it nevertheless 
failed to reproduce the data. A possible reason for this might be that Laget used in his calculations of the 7c? — > ir°np 
process the well-known Blomqvist-Laget (BL) parametrization of the pion photoproduction amplitude on the 
nucleon. This parametrization gives a good fit to the amplitude of charged pion photoproduction. But it does not 
provide a satisfactory description of it production from the proton. Because data on ir° production from the neutron 
are absent there is no possibility to check the reliability of the BL model in the description of this channel. An attempt 
to remedy this defect made in Ref. H led to a n° photoproduction operator which is not very suitable for the use in 
nuclear calculations. 

This unsatisfactory situation has stimulated a number of new theoretical investigations of the reaction d(j, ir)NN 
HHEllIlHQEEIHIll- In Ref - 01: a computation of the differential and total cross sections for this 
rocess in the first resonance region was presented. The main difference between the approaches from Refs. and 
was, that a more realistic version of an elementary pion photoproduction operator was used in the former. It 
was taken in the standard CGLN form [19( with four partial amplitudes Fi calculated with the use of the SAID 
pfij and MAID [2lJ multipole analyses. The model provided a satisfactory description of the data from Refs. 0,0)- 
Unfortunately, in Ref. H there was an error in coding the amplitude for the charged channels so that the reasonable 
description of data on the 7r _ channel should be considered as accidental. 

In a series of articles by Darwish et al. 0, ITTl IT2I Il3l Il4l ITol fl6| , the study of incoherent pion photoproduction in 
the A-resonance region was continued. ThTau?io~u^d~e~ementary operator proposed in Ref. Q which is quite 
similar to the BL operator except for slight differences in parameter values. Reasonable description of the available 
data on the total and differential cross sections was achieved in Ref. Ill[. For the first time, an attempt to analyse 
polarization observables in the reaction ^(7, tt)NN was made in Refs. |l2l Il3l IbH . ITR ITfil |. The beam asymmetry £ 
for linearly polarized photons, target asymmetries Tim, and beam-target asymmetries were discussed in those papers. 
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However, many conclusions drawn at the analysis of the polarizations were wrong as it was explained in full detail in 
Ref. H3. 

One more analysis of this process was presented in Ref. |17| where formal expressions for observables in incoherent 
pion photoproduction were given, as well as in a subsequent article 01- Using the MAID model for an elementary 
production operator, the authors studied the inclusive reaction from threshold up to 1 GeV. They obtained quite 
satisfactory agreement with the data similar to that achieved in Refs. 0, ITU, ITU . ITsl flfil | although considerable 
deviation from the predictions of the latter articles was found for many polarization asymmetries. A part of the 
deviation is a consequence of the use of wrong formal expressions in those articles. But in the cases when the right 
expressions were used, the origin of the deviation remained to be unclear. 

It should be noted that in neither of the works |n], 0, 0, ITU ITHl flit IT^ | , the question on model dependencies of 
the calculations was studied. Only the sensitivity of the predictions to the choice of the model for NN interaction was 
investigated and was found to be very small. An analogous result was also reported in Ref. 9J. But the dependencies 
mentioned do exist. They mainly stem from the elementary production operator. First, two now available multipolc 
analyses, SAID and MAID, are not equivalent and give different results for observables. The size of the deviation 
depends on the kinematic region and on the observable under consideration. For example, the total cross sections 
produced by the SAID SM04K solution at the total energy of 1232 MeV are 267 fib, 270 fib, 212 fib, and 240 fib 
for the ir°p, ir a n, ir + n, and ir~p channels, respectively. The corresponding numbers given by the MAID03 solution 
are 279 fib, 281 fib, 215 fib, and 244 fib. It is clear that this sensitivity to the choice of the analysis will be seen in 
the reaction d(-f,7r)NN too. Even different solutions for the same analysis (e.g., MAID00 or MAID03) give different 
results that leads to additional model dependence in this reaction. Second, both the SAID and MAID models are the 
parametrizations of an on-shell production operator. The latter depends on four invariant amplitudes |Tflj . Different 
options are possible for these amplitudes that are equivalent in the on-shell case. However, this equivalence is broken 
in deuteron (or more generally, in nuclear) calculations when one nucleon or both of them are off their mass shells. 
Though it is difficult to give precise numerical account for the off-shell effects, we should at least estimate possible 
uncertainties of the results introduced by them. 

We, therefore, motivate the present work by the following reasons. First, we would like to compare our results with 
those from Refs. [Til fl^L IT3I ITU ITU Il6| and Ref. in the framework of an analogous approach, i.e. in the impulse 
approximation with FSI effects, because any new calculation can serve as an independent check of which of the above 
models is valid. Second, we want to estimate possible uncertainties in theoretical predictions for observables. These 
uncertainties should be kept in mind when extracting information on the amplitude of pion photoproduction on the 
neutron from deuteron data. 

The paper is organized as follows. In Sect. ^ the kinematic relations used for the calculations and definitions 
for observables are briefly reviewed. A description of the theoretical model and its ingredients is given in sect, nni 
Section IIVI contains the results on the differential and total cross sections as well as beam and target asymmetries 
with a special emphasis on possible uncertainties of the results. In the same section we also compare our pre dictions 
with data available in the considered kinematic region and with results of other approaches given in Refs. unmm 
ITU ITHl ITfil . fl8l | . Different parametrizations for the elementary photoproduction operator are presented in Appendix lAl 



II. KINEMATICS AND DEFINITIONS OF OBSERVABLES 

Let us denote by k = (k , k), Pd = {£d,Pd), q — {s-kjO), Pi — (^liPi) an d P2 = (£2,^2) the four-momenta of 
the initial photon and deuteron and the final pion and nucleons, respectively. A symbol E 7 is reserved for the lab 
photon energy (fcP a6 = E^) and a symbol lo will be used for the photon energy in the 7c? center-of-mass (cm.) frame 
(k% n =uj = EjM/Wjd) with W yd = y/M 2 + 2ME 1 and M being the deuteron mass. 

We take as independent kinematic variables the photon energy and pion momentum q in the used frame of reference 
(generally, the lab or cm. frame) and the angles and <pp of one of the nucleons in the cm. frame of the final 
nucleon- nucleon pair. Using the equality 



W NN = 2e P = 2\j? 2 +m 2 = y/(k + Pd - q) 2 , (1) 

where m is the nucleon mass, one can find the momentum P. After boosting the momenta P and — P with the 
velocity (k+p d — q)/(k° + s d — e n ) the momenta of the outgoing nucleons are obtained and, therefore, the kinematics 
is totally determined. 

The differential cross section is given by 

da 1 m 2 £d\P\ 1 



dqdttp (2n) 5 8k ■ pd £ T £p 6 



\(m 2 m 1 \T\Xm d )\ 2 , (2) 

ni2'>ni\m ( i 



where m<z, mi, A, and m d are spin states of the two nucleons, photon, and deuteron, respectively. To obtain the 
inclusive differential cross section da/dfl^, the right-hand-side (rhs) of Eq. has to be integrated over the value of 
the pion momentum q = \q\ and the solid angle Qp 

,la fqma }d q I dn P -^- = \ s, (3) 



dVL^ J„min J dqdVL p 6 



where S is defined as 



S = J" fq 2 dq J dQp K^milTlAmrf)! 2 , (4) 



with 



f 1 ™ 2 e d \P\ () 
J (2tt) 5 8k ■ Pd e^ep ' 1 ' 

An extra factor of 1/2 must be included in the rhs of Eq. J2J in case of charged pion photoproduction. The maximum 
value q max can be found from Eq. Q at Wnn = 2m. In the cm. frame it is given by 



r ax = 2^y^ 1 (2w + ^ 2][w " d ~ (2m /i)2] ' 

q mm = 0, (6) 
where fj, is the pion mass. In the lab frame one has 

q max = (Z mQa; (e T ) = i \aE lZ + {E 1 + M)^a 2 -b^ 2 



qmin = g ™»( 07r ) = min(0, - [aE^z - {E 1 + M)^a 2 ~b^ 2 J ), (7) 

where a = (W* d - 4m 2 + ^ 2 )/2 and b = (E y + M) 2 - E 2 z 2 with z = cosG^. Note that the inequality q min ^ can 
take place only for 9^ < 90°, and at threshold energies 



4m 2 - (M - nf 
2(M - fi) 



E 1 < E™ ax = ^ (8) 



The energy E™ ax is equal to 142.6 MeV (if one takes m = (m p + m„)/2 = 938.9 MeV), 149.0 MeV, and 146.2 MeV 
for 7T°, 7r + , and 7r~ channels, respectively. Therefore, in the considered energy region, q min is equal to zero also in 
the lab frame. 

Apart from the differential cross section, single polarization observables will be considered in the article, namely, 
the photon beam asymmetry £ and target asymmetries Tim- Below, we give their definitions through the reaction 
amplitude (see also Ref. 01 )■ The photon asymmetry is 

where (do j df^)"^ is the inclusive cross section for the photons polarized parallel (perpendicular) to the zz-plane. 
Note that the minus sign in the rhs of Eq. Q is absent in the corresponding formulas from Refs. |12L ll3L ll4L UJl lm | . 



The deuteron vector asymmetry Tn and tensor asymmetries T 2 m are as follows 1 

dflp Im Yi, ((m 2 mi|T|A + l)-(m 2 mi|T|A-l))(m 2 mi|r|AO)*, 

drip (l(^2m 1 |r|A-l)| 2 + |(m 2 m 1 |T|A + l)| 2 -2|(TO 2 TO 1 |r|A0)| 2 ), 

m 2 r?i i A 

dn p Re Y ({m2mi\T\\ - 1) - (m 2 m 1 \T\X + 1)) (m 2 mi|T|A0)*, 

Re ^ (m 2 mi|T|A- l)(m 2 mi|T|A + l)*. (10) 

m^mi A 

III. THE THEORETICAL MODEL FOR INCLUSIVE PION PHOTOPRODUCTION ON THE 

DEUTERON 

The diagrammatic approach is exploited to calculate the amplitude (m 2 mi|T|Amrf) in Eq. (0). In comparison with 
Refs. 00- we reduce the set of diagrams under consideration. For example, in Ref. [g] where the threshold region was 
considered, a two loop diagram which includes simultaneously np and irN interactions had to be taken into account. 
Such a diagram is of importance at threshold energies because it involves a block with charged pion photoproduction 
from the nuclcon. With increasing photon energy this diagram becomes less important as it was shown in Ref. |6(. 
Above 200 MeV it can safely be disregarded. It is known (see Refs. 0,113) that there are kinematic regions where a 
one loop diagram with nN rescattering noticeable contributes to the amplitude. But this rather concerns the exclusive 
process 7<i — > irNN. We have checked that irN rescattering changes the final results in the first resonance region by 
only a few percentages. 

y / 7i y -~i _ / n 
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FIG. 1: Diagrams considered in this work. Two other diagrams with the permutation 1 2 are assumed. 

As a result, we retain in our calculations the two diagrams shown in Fig. ^ The pole diagram ^a) must be 
taken into account because at the integrations in Eq. 0, one goes through the kinematic regions where the relative 
momentum (pi — p 2 )/2 is small and, therefore, the deuteron wave function (DWF) has its maximum. These are the 
so-called quasifree regions. The exclusive cross section has sharp peaks in these regions. Here the inclusive cross 
section from the pole diagrams is mainly saturated. It is worth mentioning that in the peak regions the active nucleon 
Ni is almost on its mass shell. In the center of the peaks the difference between on-shell and off-shell energies of 
this nucleon is equal to the deuteron binding energy A = 2.2 MeV. Therefore, the use of the on-shell parametrization 
for a pion photoproduction operator is justified when considering the diagram in Fig. ^a). Nevertheless, as shown 
below, the off-shell dependence of calculated observables does exist even when one considers contributions from the 
pole diagram in Fig. ^a). 

Another important mechanism is displayed in Fig. IHb) . When at the mentioned integration the relative momentum 
of the outgoing nucleons, p out — (p 2 — pi)/2, decreases there are peaks in the exclusive cross sections because of 
strong final state NN interaction in the s waves (see, e.g., Refs. 0, 0, H^)- The peaks reveal themselves in a 
big contribution of diagram \2jo) to the inclusive cross section. The effect of this diagram is expected to be most 





1 The opposite sign for Tii is used in Refs. jldllfl . 



pronounced at small pion angles because in this case the low-momentum regime simultaneously for both DWF and 
NN scattering amplitude is kinematically permitted. The possibility of using the on-shell parameterization for the 
pion photoproduction amplitude at evaluation of the diagram ^b) is less evident and it is discussed below. 

Let us now write out the matrix elements corresponding to the diagrams in Fig. (see also Refs. [3. ItiL llHHsl 12*2^ ) . 
One has for the pole diagram in Fig. [Ja) 

(m 2 m 1 \T a {p 27 p 1 ,q;k)\\m d ) =^$™; Al (p 2 - yj (m 1 \T j ^ 7rNi {p 1 ,q;p 1 ,k)\Xm 1 ) 7 (11) 

where {"^(^ — Pd/2) is DWF and (mi \T J f Ii ^ TrNi (px, 9; Pi, k)\Xrhi) is the amplitude of pion photoproduction on 
the nucleon. There is one more pole diagram identical to that in Fig.^a) but with the replacement 1 <-> 2. In case of 
7T° production the corresponding matrix element should be added to Eq. 1|11[) . For the charged channels a subtraction 
of two matrix elements should be done. 



The deuteron wave function reads 



1 



(12) 



where ms = m 2 + mi, rriL = TUd — ms; Y 2 mL (p) are the spherical harmonics and Cj^f 7 fl/f are the Clebsch- 
Gordan coefficients. The present calculation is done using the CD-Bonn potential from Ref. |23] where analytical 
parametrizations of the s and d amplitudes of DWF (u(p) and w(p), respectively) is given. We note that our results 
are practically independent of the choice of a model for TV TV-interaction. Calculations with three OBEPR versions 
of the Bonn potential 0, H3 or with a separable representation [2^| of the Paris potential give almost the same 
predictions for observables. 

If it is not stated otherwise, all the results below are obtained with the production operator parametrized via the 
invariant amplitudes A{. They are defined in Eqs. (|A1|) and (|A3|) and calculated with both the SAID and MAID 
analyses. The calculation for the SAID analysis is performed in two steps. First, the CGLN amplitudes Fi in the 
cm. frame [Eq. (|A7|l ] are found making the use of electric and magnetic multipoles predicted by the analysis. We 
do not give explicit expressions for Fi through the multipoles because they are very well known. Second, one finds 
the amplitudes Ai using the relation (|A9|) between Ai and Fi. The MAID group provides users directly with the 
amplitudes Fi for the MAIDOO solution and with both Fi and Ai for the MAID03 solution. 

Strictly speaking there are other possible options for the invariant amplitudes. In particular, in Ref. another set 
of those, A' i: as defined in Eqs. ijAllfl and HA12(1 was exploited. The production operators given by the amplitudes Ai 
and A\ are equivalent in the case of on-shell nucleons as is explained in some detail in Appendix^] This equivalence 
is destroyed when the nucleons are off their mass shells. Because in deuteron calculations one deals with off-shell 
nucleons, we expect our results to be dependent on the parameterization of the elementary operator. 

The matrix element corresponding to diagram ^b) is 

, \rrbi ,m\ \ f d3 P s (Pout,rn 2 m 1 \T NN \p m ,m s m[)(m s m' 1 \T a (p s ,p' 1 ,q;k)\Xm d ) 
(m 2mi \T (p 2 ,p u q;k)\Xm d )=-m / —- ^ - 2 - 2 . (13) 

J \ Z7T ) ms f h > i Pin Pout <u 

The amplitude (m s m' 1 \T a (p Sl p' 1 ,q]k)\Xmd) in Eq. I|13[) is the same as that in Eq. I|llf) but with the replacements 
2 — > s and p\ — > p[ . The second pole diagram mentioned above with 1 <-> 2 must also be included in the integrand of 
Eq. Ijl3|l . The choice of the energy of the off-shell nucleon iVi is discussed at the end of this section. 

The half-off-shell NN scattering amplitude {p ou t, m.2mi|T/viv \Pin, i r n s m' 1 } depends on the relative off-shell momen- 
tum of the NiN 2 pair before scattering, pi„ = p s — (pi +p 2 )/2, and the relative on-shell momentum after scattering, 
Pout = (p2 -pi)/2, as 

(Pouum.m^Mp^m^) = (2*)%/%^/% Y fe^V.Cs^C^. 




x i L ~ L Y^*{p m )Y^' (p out )R J A(p out ,Pin), (14) 



where ms — m s + rh^, m' s = m 2 + mi, = mj — ms, and m^i — mj — m' s . The factors \J s ou t/m and ye^n/rn 

[sout — v Pout + m2 an( i £ in = \J Pin + 171,2 ) come from the so-called minimal relativity. The half-off-shell partial 
amplitudes R'[f L (p ou t,Pin) were obtained by solving the Lippmann-Schwinger equation for the CD-Bonn potential. 
The procedure for obtaining these amplitudes is quite direct for np and nn interactions. It should be, however, 
modified in the case of pp interaction. First, the Coulomb interaction has to be added to pure nuclear interaction. A 



method to handle Coulomb interaction in momentum space was proposed by Vincent and Phatak |27|. We do not 
discuss it here because it is described in full detail in that paper (see also Refs. [23, 13, Ull ) ■ We mention only that 
the method was applied to the 1 S'o partial wave. All other waves with J = and 1 are taken for the switched off 
Coulomb potential, ft makes no sense to include the Coulomb modifications for the waves other than 1 Sq because 
even the contribution of this latter to the observables was found to be small. As a next step, we used a prescription 
from Ref. 30] consisting of the following parameterization of the half-off-shell Sq partial amplitude for pp scattering 

Roff (Pout, Pin) = P T+»2 Ron" (Pout, Pout), (15) 
Pin ' ' 

with (3—1.2 fm _1 . The on-shell amplitude R O n {p out, Pout) is obtained with the use of the Vincent and Phatak 
method with switched on Coulomb interaction. 

All partial waves with the total angular momentum J < 3 were retained in Eq. (|14|l . In fact, however, only one 
wave, 3 Si, in the case ir° photoproduction is of importance. All other waves contribute give a few percentages to 
observables. Further details of the computations of Eq. (|13|) can be found in Ref. |5j- 
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FIG. 2: Differential cross section for the reaction d(7, n°)np in the cm. frame at 250 and 350 MeV. The dotted curves are the 
results without FSI. Only the "on-shell part" of the contribution of the diagram in Fig. 0^b) is retained in the dashed curves. 
Addition of the "off-shell part" gives the solid curves. 



Again the question emerges whether the on-shell parametrization for the pion photoproduction amplitude is appli- 
cable at calculations of the matrix element (|13|l . To discuss this point, we poin t out that at the evaluation of the 
integral in Eq. (|13JI it is assumed, in accordance with a finding from Refs. fil l22j|. the spectator nucleon N s to be on 
its mass shell. This means that the nucleon Ni momentum is —p s and its energy is equal to 2m — A — ^p 2 s + m 2 (for 
simplicity we consider the lab frame). The integral in Eq. (|I3|) is saturated at momenta p s ~ \frnK when DWF has its 
maximum. In other words, the energy of the nucleon N\ is effectively off its on-shell value by only few multiplicities 
of binding energies A. Furthermore, using the symbolic equality 

I Mr ... 1 

0(Pin -Pout) +P -Q o~ , (16) 



2 2 -n o w \rin rout) 1 J 2 2 

Pin ~~ Pout ~ zU *Pout Pi n ~ Pout 

one can split the matrix element (|13fl in its on-shell part and its off-shell part corresponding to the first and second 
terms in the rhs of Eq. l|lrij) . respectively. With only the former included we calculated the cross sections and found it 
to give the main contribution (see Fig. EJ) . Taking into account that the nucleon N s is on its mass shell, one concludes 
that this part corresponds to the case when the nucleon N[ is also on its mass shell. Therefore, the contribution of 
the amplitude Eq. I|I3|) to the inclusive cross section comes mainly from the kinematic domains in the integrand of 
Eq. (|13JI where the nucleons Ni and N[ are close to their mass shells. The same conclusion holds true also for other 
observables. Therefore, the on-shell parameterization for an elementary pion photoproduction operator is applicable 
in this integrand. Nevertheless, some dependence of the FSI amplitude [Eq. i|13|) ] on off-shell effects is expected. 

As in Refs. 0, @, all summations over polarizations of the particles in Eqs. (|ll|l and 1I3|) as well as the three- 
dimensional integration in Eq. (|13|l have been carried out numerically. The number of chosen nodes at this integration 
and that in Eq. @ was taken to be sufficient for prediction of observables with the numerical accuracy better than 
2%. 



IV. RESULTS AND DISCUSSION 



A. Differential and total cross sections 



We begin our discussion with the results for the neutral channel. In Fig. the predicted differential cross sections 
of 7T° production are shown in the energy region between 208 and 419 MeV together with experimental results from 
Refs. 0,0 2 ■ The displayed cross sections are obtained with a pion photoproduction operator parameterized through 
the amplitudes At and the SAID FA04K solution. One can see one more confirmation of a prediction from Refs. 0>l3 
that the effect of np final state interaction should lead to a reduction of the cross section and this reduction is the 
stronger the smaller the pion angles are. This effect is mainly attributed to the strong np interaction in the 3 «Si wave. 
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FIG. 3: Differential cross section for 7r° production in the photon-nucleon cm. frame obtained with the FA04K solution of the 
SAID analysis and with the parametrization of the operator through the invariant amplitudes Ai. The dotted (solid) curves 
are our predictions without (with) FSI. Data are from Ref. (•) and Ref. (o). 



Without FSI the model completely fails to reproduce the data. After including FSI, the curves move to the data 
points although a reasonable description of them still remains to be achieved. At 60° < 9;^ < 120°, the predicted 
cross sections overestimate the data by about 10 — 20%. 

Possible explanations of the disagreement between the data and the present model can be looked for in the elemen- 
tary photoproduction operator. As is explained above, the on-shell parameterization for the latter can be used when 
studying the inclusive channels. Nevertheless, different representations of the operator, which are equivalent in the 
on-shell case, turn out to be not quite equivalent when one or two nucleons are off their mass shells. Because the so- 
phisticated phenomenological analyses like MAID or SAID provide the elementary amplitude for the on-shell nucleons 
only, it is hardly possible to give precise quantitative account for the off-shell effects using these amplitudes. One can, 



2 In Refs. the differential cross sections are given in the so-called "photon-nucleon cm. frame" . R elations needed to transform the 

cross sections and angles from the -yd cm. frame to the frame mentioned are presented in Refs. Isl llOl . 



however, estimate the possible size of the off-shell effects by performing calculations with different representations of 
the operator. As an example, we make use of two forms of the operator given in Appendix corresponding to the 
amplitudes Ai (|A1|) and A\ I|A11(I . In addition, the MAID03 solution was used to parametrize this operator. 
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FIG. 4: Differential cross section for the 7r° channel in the photon-nucleon cm. frame with different parametrizations of the 
production operator: dashed and solid curves are obtained with the amplitudes A' t and Ai, respectively, for the SAID FA04K 
solution. Dotted curves are obtained with the amplitudes Ai and the MAID03 solution. Data described in the legend to Fig. [3] 



The cross sections, shown in Fig. exhibit quite noticeable dependencies to the different parametrizations of the 
operator. At 208 MeV, the cross section is sensitive both to the choice of the analysis and the on-shell form of the 
operator. With increasing photon energy, the sensitivity to the analysis diminishes but the sensitivity to the form 
of the operator remains to be noticeable. The regions overlapped by the curves can be considered as characterizing 
the size of possible uncertainties introduced by the pion photoproduction operator. They should be kept in mind at 
attempts to extract the cross sections on the neutron from deuteron data. 
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FIG. 5: Differential cross section for the reaction d(-y,-K°)np in the photon-nucleon cm. frame obtained in three models: the 
present calculation with the amplitudes A[ and the MAID03 solution (solid), Ref. [TH (dotted), and Ref. |isj] (dashed). Data 
described in the legend to Fig. 

A comparison of our results with those from recent works [ll| and [18j is presented in Fig. Because in the 
latter article a photoproduction operator was parametrized via the amplitudes A\ and the MAID03 solutiom we 
give the comparison with the same operator. One can see a satisfactory agreement with the results from Ref. 18]. 
Slight deviation might be attributed to the use of the different parametrizations for the half-off-shell NN scattering 
amplitude. At the same time, our differential cross section exhibits quite a different behavior compared to that from 
Ref. • Reasons responsible for this incongruity of the results can be in the use of both different half-off-shell AW 
scattering amplitudes and elementary photoproduction operators. We suppose the latter reason to be more probable. 
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FIG. 6: Total cross section for the reaction d(y, Tv°)np. (Left panel) Results with the FA04K solution of the SAID analysis 
and with the parametrization of a pion production operator through the invariant amplitudes Ai. The dotted (solid) curves 
are the predictions without (with) FSI. (Middle panel) Dashed and solid curves are obtained with the amplitudes A\ and Ai, 
respectively, for the SAID FA04K solution. The dotted curve is obtained with the amplitudes Ai and the MAID03 solution. 
(Right panel) The present calculation with the amplitudes A\ and the MAID03 solution (solid), Ref. [Till (dotted), and Ref. [Ts| 
(dashed). Data described in the legend to Fig. 
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FIG. 7: Differential cross section for the reaction d("/, n~)pp in the lab frame obtained with the amplitudes Ai and the SAID 
FA04K solution. The dotted curves are contributions from one of the pole diagrams in Fig. Successive addition of the second 
pole diagram and FSI leads to dashed and solid curves, respectively. Data described in the legend toRef. [3ll |. 



After integrating Eq. J2J) over the solid pion angle one obtains the total cross section for a given channel. In Fig. 
the total cross section for 7r° photoproduction is shown. It is clear that everything told above on the differential 
cross section also holds true for the total cross section. In particular, one can see that the model without FSI clearly 
overestimates the data. Inclusion of FSI strongly reduces the cross section although it still noticeable overestimate 
the data from Refs. and for all parametrizations of an elementary photoproduction operator. Our total cross 
sections are about 5% higher as compared to those from Ref. in the vicinity of the peak. Note also that in 
comparison to the prediction from Ref. [Tl|, our peak position is shifted on about 10 MeV to smaller energies that 



seemingly is because different elementary photoproduction operators used in two models. 
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FIG. 8: Differential cross section for the reaction d('y,n~)pp in the lab frame with different parametrizations of a production 
operator. Notations of the curves as in Fig. |I] Data are from Ref. [3l| . 



In discussion of charged pion photoproduction we restrict ourselves to the d("f, ir~)pp channel only because results 
for the (i(7,7r + )nn channel are very similar. The dotted curves in Fig. that correspond to the contribution of one 
pole diagram, reproduce the behavior of the angular dependence for the differential cross section of the elementary 
reaction Tt~p. In particular, at energies above 400 MeV, a sharp peak at forward angles because of the pion 

exchange in the t channel is clearly seen. One can see that at 0„- > 90° the cross section from two pole diagrams 
is practically equal to twice the cross section from one diagram. The reason for this is that at backward angles the 
events where both nucleons have small momenta correspond to the high momentum components of the deuteron 
wave function and, therefore, both diagrams cannot work in the quasifree regime at the same time. As a result, 
the interference term is very small for backward angles. Of course, this conclusion is valid for all channels. In full 
agreement with a finding from Refs. [IlllTsj] . the effect from FSI has only a marginal impact on the differential cross 
section. 

As is seen in Fig. |Hl the sensitivity of the cross section to the different parametrizations of the photoproduction 
operator is not as strong as in the case of the 7r° channel. It remains to be visible only at forward angles. 
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FIG. 9: Differential cross section for the reaction d(7, ir~)pp in the lab frame obtained in three models: the present calculation 
with the amplitudes A[ and the MAID03 solution (solid), Ref. (dotted), and Ref. [Isj (dashed). Data are from Ref. |3lj |. 

A comparison of our results with those from recent works [Til Il8| is presented in Fig. [5] One can see the good 
agreement with the results from Ref. ,18.]. We, however, have expected better agreement because the FSI effect is small 
for the charged channels and, as stated in Sect. IIIII the results are independent of a choice of DWF. The deviation 
at forward angles is even somewhat bigger because the Coulomb forces between the protons have been disregarded in 
Ref. [is|. Their effect consists in the decrease of the cross section on 5 — 10% at zero angle in the energy region from 
250 to 500 MeV and becomes to be negligible at 0,r > 30°. Therefore, reasons for the slight deviations between our 
results and these from Ref. 0] remain to be inve stig ated. The disagreement with the results from Ref. ^lj is also 
difficult to explain because, as it is stated in Ref. [llj, the operator used in that work provides the differential cross 
sections of the elementary reaction on the nucleon close to those given by the MAID analysis. 

The total cross section for 7r _ photoproduction is shown in Fig. 1101 Here the FSI contribution is much smaller than 
that for 7T° production and leads to a slight decrease of the cross section above 350 MeV. We find satisfactory agreement 
with data from Refs. [3ll l32L l33|. At the same time a data point from Ref. [34| at 250 MeV lies markedly below 
both our predictions and data from Refs. 0>IH2- The sensitivity of the results to the choice of the photoproduction 



operator is rather small. A comparison to the results from Refs. [TH Il8[ shows that three models give very similar 
results. 
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FIG. 10: Total cross section for the reaction d(7, n~)pp. Notations of the curves as in Fig. E] Data are from Refs. [3lj (■), 
P (A), HI (□), and (.). 



Having results for the total cross sections in all the channels mentioned above one can try to make predictions 
for the total photoabsorption cross section on the deuteron in the first resonance region. Of course, two more 
reactions contribute to it as well. These are coherent ir° photoproduction from the deuteron, 7c? — > n°d, and deuteron 
photodisintegration, 7c? — ► np. Predictions for the former are taken from a model built in Ref. j3^ which provides a 
good description of data from Ref. pj. The total cross section for the latter reaction is calculated making use of a 
phenomenological fit |3j| to available experimental data on deuteron photodisintegration up to 440 MeV. 
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FIG. 11: Total photoabsorption cross section per nucleon for the deuteron from 150 to 500 MeV. Contribution of the reaction 
jd — > Tv°d is shown as dotted curve. Contribution from deuteron photodisintegration is included in dashed curve. Filled area 
includes contributions from all channels 7c! — > irNN (see text). Data are from Refs. |37| (o) and |3S| (•). 

In Fig. 1111 we present our results for the total photoabsorption cross section per nucleon for the deuteron. The 
filled area includes the uncertainties discussed above, because os the variations of the elementary photoproduction 
operator. It is seen that the predictions even with allowance for these uncertainties are noticeable above the data from 
Refs. [37], in the peak region. In the center of the peak at about 320 MeV we find our prediction of (543 ± 7) /ib 
strongly overestimating the experimental value of (452 ± 5) fib. We have no explanation for this disagreement. There 
are reasons for the assumption that the total cross sections from Refs. [37i l33| may be too low. Such an assumption 
is supported by the study of other electromagnetic processes. As an example, one can mention that the sum of 
the electric and magnetic polarizabilities of the neutron calculated with those cross sections turns out to be notably 
underestimated (see Ref. |3|| for a more detailed discussion). 




B. Beam asymmetry for linearly polarized photons 

The beam asymmetry £ for 7r° production at three selected energies of 250, 350, and 500 MeV is displayed in 
Fig. E| In IA it is negative at all energies. FSI results in a decrease of the magnitude of S. The influence from 
FSI is noticeable at the lowest energy and forward angles. With increasing energy the effect of FSI becomes smaller 
although not negligible even at the highest energy. 

As in case of the differential cross section, the beam asymmetry depends strongly on the form of the elementary 
production operator. As is seen in Fig. 1121 at 250 MeV there exists noticeable sensitivity of S to the choice of the 
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FIG. 12: Angular distribution of the asymmetry E for the reaction d(j, ir°)np. (Upper panel) The dotted (solid) curves are our 
predictions without (with) FSI obtained with the SAID FA04K solution and with the amplitudes Ai. (Middle panel) Dashed 
and solid curves are obtained with the amplitudes A\ and Ai, respectively, for the SAID FA04K solution. Dotted curves are 
obtained with the amplitudes Ai and the MAID03 solution. (Lower panel) The present calculation with the amplitudes A\ and 
the MAID03 solution (solid), Refs. (dotted), and Ref. EH (dashed). 






analysis of photomeson amplitudes and. to the choice of their representations in terms of Ai or A[ . This sensitivity 
mainly reflects the difference between the parametrizations FA04K and MAID03 in case of 7r° production near 250 
MeV. In the A region and at higher energies, the asymmetry is practically independent of the analysis. But here 
the influence of the choice of the amplitudes Ai or A[ becomes visible. Only at 350 MeV these amplitudes lead to 
very close results for the asymmetry. From this finding one may conclude that the A region is promising for a model 
independent determination of E for the ir°n channel from deuteron data. Outside this region, results of such an 
extraction will be strongly model dependent. 

A comparison of our results with those from Refs. 0, 0] and E3l * s a l so presented in Fig. Our predictions 
are similar to those from the latter work although there is a disagreement at forward angles for 250 MeV. As to the 
predictions from Ref. [l2l ITsT l , we have a notable disagreement with the results from those works in absolute size of E 
at 250 and 350 MeV and in the form of the angular distribution at 500 MeV. One should mention a strange result of 
Refs. [l2,Ilil consisting in the statement that E does not vanish at = and it, as it has to be because of helicity 
conservation [r^.ll8| . 
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FIG. 13: Angular distribution of the asymmetry E for the reaction d(7,7r°)np. Notation of the curves as in the middle panel 
of Fie. IT21 Data (preliminary) are from Ref. 40]. 



Recently first preliminary data on the asymmetry E in the tt° and it channels at a few energies between 265 and 330 



MeV have been reported by the LEGS collaboration [40( . A comparison with the data for ir° production is presented 
in Fig. 1131 One can readily see a satisfactory agreement with the experimental values for all parametrizations of the 
photoproduction operator. Only near 90° the asymmetry is slightly overestimated in absolute size when using the 
operator built with the SAID model. 
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FIG. 14: Angular distribution of the asymmetry £ for the reaction d(-y,n~)pp. Notation of the curves as in Fig.[H Results 
from Refs. [TflL H3 | are shown as dotted curves of the lower panel. 



The beam asymmetry E for ir~ production is shown in Fig. 1141 Remaining negative as in the case of tt° production, 
it shows a quite different behavior. One can see a sharp peak near 9 T ~ 30°. The effect from FSI is much smaller 
and has a noticeable impact on £ only at the lowest energy. At the highest energy it is negligible. The SAID and 
MAID analyses give quite different results for £ at the lowest energy. This difference is diminishing when the energy 
increases. Notable influence on the predictions at 250 MeV has also the form of the production operator. 

In the same figure we compare our predictions to those from Refs. Note that polarization observables for 

the 7T~ channel were calculated also in Refs. [l2l[l^ |. In some cases they are in notable disagreement with those given 
in Refs. |l5lll6| |. Authors do not explain reasons for the deviation. In this situation we preferred to make a comparison 
with predictions from the more recent works. Our results are in reasonable agreement with the predictions of Ref. |18| . 
but in substantial disagreement with the results of Refs. ^J^l both for the form of the angular distribution and for 
the absolute size of £ 3 . The disagreement is drastic at the highest energy and can hardly be caused by the use of a 
different elementary production operator. It is likely, that there is an unnoticed computational error in Refs. |l5lll6|. 
leading also to the odd results of nonzero asymmetry at = and n. 



3 In Refs. |3 E3 the angular distribution of S is given at E 1 = 200, 270, 330, 370, 420, and 500 MeV. The curves at 250 and 350 MeV 
displayed in Fig. 1141 have been obtained by a quadratic interpolation. The same procedure is used to obtain the target asymmetries 
Tim f°r tt — production discussed below. 
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FIG. 15: Angular distribution of the asymm etry E for the reaction d(7, tv )pp. Notation of the curves as in the middle panel 
of Fig. [TTI Data (preliminary) are from Ref. |40|. 



In Fig. E| we compare our results with the preliminary data from the LEGS collaboration at 270 and 330 MeV. 
Although some model dependence of the predictions exists, the calculated asymmetry is too small in its absolute size 
at all the parametrizations of the production operator. 

We do not consider here the tt + channel because all the conclusions just drawn for tt~ production remain to be 
valid for ir + production as well (see also Ref. [l8j) and move to a discussion of the target asymmetries. 



C. Target asymmetries for polarized deuterons 

The target asymmetries Tim for ir° and ir~ production are shown in Figs. 1161231 The asymmetries for ir + production 
are very similar to those in the 7r~ channel and are not discussed below. 

In accordance with results from Ref. we have found that the form the angular distribution of Tyi in the 7r° 
channel changes notable with increasing energy. FSI effects are rather small. The sensitivity of Tu to the choice 
of the photoproduction operator is quite small at the lowest energy. This is not the case at higher energies. The 
vector asymmetry is sensitive to both the choice of the analysis and the form of the operator. A good agreement of 
our results with those from Ref. 0] is seen. Only at forward angles < 30° we have found some deviation. The 
predictions from Refs. (l2L Il3| totally contradict our results both in the form of the angular distribution and in the 
absolute size at 250 and 350 MeV 4 . For instance, we observe a maximum around 130° — 140° at 250 and 350 MeV 
with Tn ~ 0.4 in the center of the peak, whereas in Refs. [T2I H"^ the vector asymmetry is close to zero at these 
energies for Q w > 90°. We have no a reasonable explanation for this disagreement. Only at the highest energy one 
has a good agreement between the three calculations. 

Figure lTTI shows that FSI effects on the vector asymmetry for ir~ production are much smaller than for 7r° production. 
An analogous result has been reported in Ref. |lj| . As in the case of the tt° channel, the vector asymmetry is practically 
independent of the photoproduction operator at the lowest energy. At higher energies we observe some sensitivity of 
Tn to that operator especially in the angular region between 30° and 120°. One can observe very good agreement 
with the results of Ref. 0] but significant disagreement with those of Ref. is evident. For instance, we do not 
find a peak near 30° at 250 MeV predicted in Ref. In the peak position our value for T\\ is by a factor of ~ 3 
smaller than that in Ref. |l6[ . 

The tensor asymmetry T20 for 7r° production is displayed in Fig. 1181 It is small in absolute size both in IA and 
IA+FSI. The FSI effect is very large at forward angles as it was previously found also in the case of other observables 
in the neutral channel. The SAID and MAID models give very close predictions for T20 but the results are sensitive 
to a form of a production operator. The asymmetry T20 is the first observable for which we have found a substantial 
deviation from results of Ref. |l8j|. Our resulting as ymm etry shows a sharp minimum at forward angles, whereas a 
sharp maximum was found in Ref. [Tsj . The model |lfij | predicts even deeper minimum at forward angles than that 
in our calculation. 

As is seen in Fig. 1191 the asymmetry T20 for the ir~ channel is forward peaked in IA. FSI has only a marginal impact 
on T20. The resulting asymmetry is very small for ©„■ > 30° especially at high energies. In this kinematic region 
the results are practically independent of the choice of the production operator. The three approaches predict close 



4 When making comparisons to results from Refs. Il2lll il7tl one should keep in mind that in those works the asymmetry T11 is defined 
with the opposite sign. 



peak values of T20. Some deviation takes place in the regions where the asymmetry is small. Note that recently the 
tensor target asymmetries in n~ photoproduction on the deuteron have been studied in Ref. [4l| . Because the authors 
considered the exclusive reaction, a comparison of our predictions with the results from that work is impossible. 

The target asymmetry T 2 i for 7r° production in Fig. |2*U1 shows drastic FSI influence. The resulting absolute size of 
T21 is small. It does not exceed 0.1 for O^ > 30° at all energies and shows a notable model dependence. There exist 
quite significant differences between the results of our model and those from Refs. and [l^- As is seen in Fig. 1211 
practically all the above conclusions remain to be valid for tt~ production too. Only at 250 MeV we observe a smaller 
influence of FSI and a good agreement with the results of Ref. is evident at this energy. 

Predictions for the target asymmetry T22 for ir° production are shown in Fig. 122 This asymmetry is very small in 
IA. Its absolute size is less than 0.03 in the kinematic region under consideration. FSI manifests itself in a pronounced 
peak around 20° although even in the center of the peak the asymmetry T22 is still small being less than 0.12. The 
sensitivity of T22 to the production operator varies with the kinematics. It is quite small at 250 and 350 MeV but 
becomes to be notable at 500 MeV. It is seen in Fig. 1221 that the differences to the results of the three approaches are 
quite significant. 

More pronounced peaks around 20° but in IA are seen in Fig. 123 where the target asymmetry T22 for 7r~ production 
is displayed. Influence of FSI is quite small and seen only at the highest energy for > 60°. One notes a sizable 
dependence of the predictions to the production operator at the lowest e nerg y. It is much smaller at higher energies. 
We find good agreement with Ref. [lj| but significant differences to Ref. [l^ are evident at 250 and 350 MeV. 
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FIG. 16: Target asymmetry Tn for the reaction d{-y, -K°)np. Notation described in the legend to Fig. 1121 



V. SUMMARY AND OUTLOOK 

In the present work we have studied incoherent pion photoproduction on the deuteron in the first resonance region 
taking into account diagrams with a plane-wave final state and with NN interaction in the final state. Particular 
emphasis has been laid on the discussion of possible uncertainties introduced into the model by the elementary 
operator of pion photoproduction on the nucleon. We have demonstrated that the use of different forms of the on-shell 
operator has a notable impact on predictions both for the unpolarized cross section and for polarization observables, 
in particular, the beam and target asymmetries. It is evident that these uncertainties will manifest themselves in the 
corresponding variations of the amplitude of pion production on the neutron extracted from deuteron data. We have 
not studied beam-target asymmetries in the present work leaving their consideration for a subsequent publication. 
However, we are confident that analogous uncertainties will be also seen in these asymmetries. 

We have also carried out a detailed comparison of our predictions for unpolarized cross sections and beam and target 
asymmetries with recent results from Refs. [121 i!3L Il5l Il6| and Ref. ^g. For most observables we have found good 
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FIG. 17: Target asymmetry Tn for the reaction d( / y,n~)pp. Notation of the curves as in Fig. 1141 only in the lower panel, 
results from Ref. Hal are shown in dotted curves. 
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FIG. 18: Target asymmetry T20 for the reaction d(7, ir°)np. Notation described in the legend to Fig. 1121 



agreement with results of the latter work. However, as a rule our predictions are in significant deviation from those 
in the former works. Of course, part of the disagreement can stem from the use of a different elementary production 
operator. But in many cases the deviation is too big to have such an explanation so that the reasons for the deviation 
are still to be understood. 

Practically for all paramctrizations of the photoproduction operator, our predictions for the unpolarized differential 
and total cross sections in the 7r° channel are too big in comparison to the available data. The agreement with data in 
the tt~ channel is quite good. The situation with the description of the preliminary data from the LEGS collaboration 
on the beam asymmetry £ is opposite. The agreement is satisfactory for ir° production but our predictions for the 
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FIG. 19: Target asymmetry T20 for the reaction d{ r y,-K~)pp. Notation described in the legend to Fig, 1171 
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FIG. 20: Target asymmetry T21 for the reaction d(7, -K°)np. Notation described in the legend to Fig. 1121 

tt~ channel clearly underestimate the experimental values in absolute size. 

An imp ortant problem to be solved is the significant disagreement between the theory and experimental data from 
Refs. |37ll38l | on the photoabsorption cross section for deuterium in the vicinity of the A peak. Even with the allowance 
for all the uncertainties considered in the present model, we have found our predictions to overestimate significantly 
the measured values. We suppose that the problem might be in the data themselves. Using the deuteron values, 
the authors of Ref. extracted the total photoabsorption cross section on the neutron that, near the A resonance 
energy, is by ~ 20% lower than that presently predicted by the SAID and MAID analyses. New measurements of the 
deuteron cross section in the peak region would be of very importance to clarify the situation. 

Much additional work should be also done to improve the theoretical model. In particular, one should try to 



E =250 MeV 




E =500 MeV 



0.0 



30 60 90 120 150 180 30 60 90 120 150 180 30 60 90 120 150 180 
0.4, , , , , 1 1 0.2, , , , , . 1 0.2 r 




30 60 90 120 150 180 30 60 90 120 150 180 30 60 90 120 150 180 
®T [cleg] [deg] 6* [deg] 

FIG. 21: Target asymmetry T21 for the reaction d('y,ir~)pp. Notation described in the legend to Fig. 1171 
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FIG. 22: Target asymmetry T22 for the reaction d(7, -K°)np. Notation described in the legend to Fig. 1121 



take into account two-loop diagrams. As has already been mentioned in Refs. |3|, |42( for the case of exclusive pion 
production, there are the kinematic regions where two-loop diagram with NN rescattering in the intermediate state 
can be of importance. However, it is not an easy task to include that diagram in the model for inclusive production 
because numerical calculations become to be extremely time consuming. Also, in view of the notable sensitivity of 
the predictions for many observables to the choice of the elementary operator, there is urgent need to develop realistic 
models which take into account of off-shell effects in pion photoproduction on bound nucleons. 
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APPENDIX A: AN ELEMENTARY PION PHOTOPRODUCTION OPERATOR 

The invariant pion photoproduction amplitude can be written as 

4 

T~/N—>irN' = u(p')^2A i (s 1 u,t)r i u(p), (Al) 

i=l 

where 

s = (k+p) 2 = {q + p'f, u={k-p'f = {q-pf, t=(k-q) 2 = (p-p') 2 , (A2) 

and 

r 1 = 175^, 

T2 = «75 [q ■ e (P + p') ■ k - q ■ k (p + p') ■ e] , 
T3 = «75 (<? ■ k(- q- af) , 

T 4 = e^j^ePk' . (A3) 
The matrix 75 and antisymmetric tensor €^ upa are fixed according to the conditions 

75 = + (io) and 60123 = +L ( A4 ) 

In the spinor form, the matrix T reads 

(m 2 |T|Ami) = (m 2 \L + ia ■ K\\ mi ). (A5) 



Contributions of the amplitudes Ai to the matrix T l|A5(l in an arbitrary frame are as follows 
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(A6) 



where S = k x e, G = f x f , E± = E ±m, E'± = E' ± m, N = x /E+/2m, and N' = x/E' + /2r 
The photoproduction operator in the cm. frame has the well-known form [l9l | 
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where W = y/s and the superscript asterisk is used for the corresponding quantities in the 7-/V cm. frame. A 
comparison of Eq. (|A6(I in the cm. frame with Eq. (|A7(I gives the following relation between the amplitudes Fi and 
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(A8) 



where W± = W ± to. The inverse relation is as follows 
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/ 



Note that the corresponding formula in Ref. |4J| contains missprints for the amplitudes A\ and A 3 . One should 
emphasize that all variables in Eqs. I|A8J) and IjAQjl are taken in the cm. frame, in particular 
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Another set of the invariant amplitudes A\ was used in Ref. 0] 
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(All) 



where 

ri = £75 = Ti, 

r' 2 = 2i7 5 (e • p'k -p- e ■ pk ■ p') = -T 2 , 

T' 3 = ij 5 (#e-p- (k-p), 

T 4 = H50e-p' -ik-p'), (A12) 

There are the following relations between the amplitudes Ai and A\ 

A\ =A X - 2mA 4 , A' 2 = -A 2 , A' 3 = -A 3 - A 4 , A' A = A 3 - A 4 . (A13) 

Using the spinor form of the amplitude T [Eq. l|A6(l ] and that for T" given in Ref. one can show that the 
difference between T and T' reads 

NN' - 

T-T'= -M ——(m 2 \6p' ■ S - S'p- S + icf- [e(5(muj + p'-k) + 5'(mLu+p- k)) + k{Sp' ■ e + S'p- e)]\Xmi), (A14) 
E + E + 

where <5 = p 2 — m 2 and 5' = p' 2 — m? . Therefore, two representations of the pion production amplitude, Eqs. (|A1|) 
and I|A11J1 . are equivalent only in the case of on-shell nucleons when S = S' = 0. 
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